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Abstract 

We present an alternative approach for calculating helicity am- 
plitudes for processes involving both massless and massive fermions. 
With this method one can easily obtain covariant expressions for the 
helicity amplitudes. The final expressions involve only four vector 
products and are independent of the basis for gamma matrices or spe- 
cific form of the spinors. We use the method to obtain the helicity 
amplitudes for several processes involving top quark production. 
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1 Introduction 



The conventional approach used to determine decay rates and cross sections 
is to square the Feynman amphtude, average over initial, and sum over final 
spin states. This approach is very inconvenient to implement if the number 
of Feynman diagrams and the number of final state particles is large. This 
is true even with the use of all the modern tools for symbohc computation. 
The conventional approach also has the drawback that in summing over spin 
states one loses information on spin correlations. Such correlations may be 
important for exploring new physics at future colliders. 

An alternative approach, which has gained popularity in the past decade, 
is to compute the helicity amplitudes numerically. By doing this, spin cor- 
relations are easily explored and the process of squaring the amplitude is 
trivial. The idea of calculating helicity amplitudes is probably as old as the 
conventional approach and might even pre-date it[l, 2]. More recently, mo- 
tivated by the need to consider more complicated processes, a number of 
novel methods[3, 4] for obtaining helicity amphtudes have been developed. 
Most of these, however, are viable only in the limit where fermion masses 
can be ignored. In this limit helicity and chirality are equivalent, leading to 
enormous simplifications; this is particularly evident in processes involving 
the VT-gauge bosons. Most of the recent so-called hehcity methods[3, 4, 5] 
take advantage of these simplifications [6]. 

Although, there have been several attempts to generalize the helicity methods 
for massless fermions to the massive case, the resulting procedure is rather 
cumbersome and requires the introduction of arbitrary reference momenta[7] . 
It is the purpose of this paper to present a simple and systematic method 
for calculating helicity amplitudes which works for processes involving both 
massless and massive fermions. The approach that we present leads to a 
covariant expression for the helicity amplitude in terms of a minimum number 
of four vector products and is independent of the basis for Dirac matrices or 
specific forms of the spinors. Thus, we believe that our approach leads to an 
analytic expression which is relatively simple and easy to compare with. 

For illustration we apply our method to compute the helicity amplitudes for 
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various processes involving the top quark. We consider the processes: 

e+e~ tt 
77 — > tt 

2 The Method 

When evaluating a Feynman amplitude involving a pair fermions and several 
gauge bosons, in the initial or final state, the amplitude is expressed as a 
sum of terms which have the form, 

u{pA)$SS2---$2n+l<P'^>!)^ (1) 

where (p. A) and (p'. A') are the momenta and helicity of the external fermions. 
The particular labeling of the indices is chosen for convenience in the ma- 
nipulations that follow. Note that, as written, (1) has an even number of 
7-matrices sandwiched between spinors. The case where the number of 7- 
matrices is odd is treated by setting = /. 

In the usual approach one squares these terms and sums over the fermion 
spins to obtain a sum over terms of the form, 

]■ (2) 

The traces can then be evaluated and the result will be a function of Minkowski 
products of the external particle four momenta. The expression is then easily 
integrated over phase space by use of Monte Carlo methods. In this approach 
one does not need explicit values of components of free spinors (which de- 
pend on the particular representation used for the 7-matrices). This method 
is also convenient because, in addition to getting rid of the 7-matrices and 
expressing the amplitude in terms of simple "dot" products, the final result 
is independent of the unobservable fermion spins. 

In complicated cases where the amplitude involves the sum over a large num- 
ber of Feynman diagrams, and each diagram involves many fermion-gauge 
boson vertices, the number of traces and the length of the traces may render 
the standard approach impractical. The standard approach is also inconve- 
nient in cases where interference effects are of interest and one may want 
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to explicitly display the relative phases between the different Feynman dia- 
grams. In these cases it is more convenient to obtain an expression for the 
amplitude corresponding to each graph which is simple to evaluate numeri- 
cally. The square of the amplitude is then obtained by just squaring a single 
complex number. 

For processes that only involve massless or nearly massless fermions there 
have been many papers where the so called helicity methods for evaluating 
tree level amphtudes are developed[3, 4, 5]. In contrast to the standard 
approach, the amplitude in the Helicity Method is expressed in terms of 
spinor products of the form 

^u(p, +)u(p',-). 

In this definition the p and p' are hght-like four-vectors. The ± refer to 
the hehcity of the fermion. Note that \s{p,p')\ ^ = p-p', so that up to a 
phase s{p,p') contains essentially the same information as p-p'. These meth- 
ods prove relatively simple to implement and lead to compact expressions 
valid whenever the fermion masses can be ignored. Generalizations of this 
approach to the massive fermion case[7, 8] exist, but the results are cum- 
bersome: the method requires the introduction of extra light-like momenta 
and the expression of the spinor for a massive particle as a sum over spinors 
for massless particles. Consequently, as the number of external particles in- 
creases, the number of terms for each diagram grows and results quickly be- 
come unmanageable and more prone to human error. In another approach [9] 
some authors found it more convenient to obtain the amplitudes by perform- 
ing all the 7-matrix products numerically. Although this approach appears 
to be efficient, the result is not expressible in a simple compact analytic form 
which one can study. We finally note that all the approaches quoted above 
require writing down particular expressions for the spinors. 

We will see below that with the approach presented here, it is possible to 
actually write the amplitude in terms of "dot" products without the use of 
particular expressions for the spinors. The procedure is relatively simple and 
leads to compact analytic expressions which are easy to compare with and 
put into code to evaluate numerically or symbolically. 

The basic idea is to observe that (1) can be re- written as the trace, 

tr[uip', A') ^ u{p, XW^i^2 ■ ■ ■^2n+i\- (3) 
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Although this trivial observation has been used by several authors in the 
past[10, 11]. there was no attempt was made to make use of the Michael- 
Bouchiat identity (as we do below) to obtain an expression for u{p', X')u{p, A). 
The object, u{p', A') ®u{p, A), is a four by four matrix in Dirac space and can 
be expressed in terms of an orthogonal basis of the four-dimensional Dirac 
space [12]. We do not follow this approach here because we found it simpler 
to make use of a spinor identity due to Michel and Bouchiat [13], 

A')«(p, X) = \{i> + m){5xx' + ld<y\x')- (4) 

The cr's are the three Pauli matrices and there is an implied sum over i — 
1,2,3. Note that in the case where A' = A expression (4) reduces to the usual 
projection operator for a state of momentum p and helicity A. The jy's are 
defined such that, 

Vi-P = 0. (5) 

They are the spin vectors corresponding to the four momentum p. For, 

p=\p\ (sin 9 cos 4>, sin 9 sin 0, cos 9) , 

a standard representation of the 77's is given by, 

rji = (0; cos^cos^, cos^sin^, — sin^) 
V2 = (Oj ~ sin (p, cos (p, 0) 

nS = (^;^^). (6) 

m m \p\ 

where rj^ = (0;5|.) for the case when \p\ = 0. 

For completeness we have included a proof of identity (4) in the appendix. 
The identity is only valid if both spinors correspond to the same momenta. To 
make use of it we must generalize it to obtain an expression for u{p', X')u(j>, A). 
li m' — m this could be done by applying a boost operator to both sides of 
(4). It is, however, simpler and more general to observe that any spinor 
u{p', A') can be expressed in the form 

u{p\\') =XJ2Cx'xi]^' + m')u{p,\), (7) 
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where //^ = m'^, and jV is a normalization factor. This expression is vahd 
whether or not m' = m. Note that the heUcity index A is defined with 
respect to the direction of p , while A' is defined with respect to p' ; i.e. the 
frames in which A and A' are defined are rotated with respect to each other. 
Multiplying (7) on the right by u{p, A) and using identity (4) we obtain, 



(8) 



with, 
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'2{p-p' + mm') 

It is easy to verify that the matrix Cyx is unitary in helicity space as it 
should be. Using (7) in (4) we obtain the desired generalization, 

uip', A') uip, A) = ^ ^ Cx'x" + m')(^ + m){5xx" + id^^lx")- 

^ A" (9) 

Since Cxx: is unitary it drops out of our expressions upon summing the square 
of the helicity amplitudes over all helicities. Thus, if in the end one is going 
to sum over all helicities one can simplify (9) by setting Cx'x" = ^x'X"- The 
resulting amplitudes will not correspond to the correct helicity amplitudes 
but the process of squaring them and summing over all helicity indices will 
indeed generate the correct result for the sum of the squares. 



If spin correlations are of interest then one needs to make use of the specific 
expressions for the Cxx'- The expressions follow from (8) 



EE'— \p\ \p'\ +mm' 
p-p' + mm' 



1/2 



e 9' 

cos - COS — 

2 2 



d 9' -111 i\ 
sm - sm — eH"^ -'^> 
2 2 



EE'+ \p\ \p'\ 



mm 



p-p' 



mm 



1/2 



e 6' 

cos - sm — e 
2 2 



. e 6' 

— sm - cos — e 
2 2 



C 



++ 



= -c 



+- 



(10) 



The angles (9, (f)) and {$', (j)') correspond to the axial and azimuthal angles 
for p and p ' respectively. In the center of momentum frame where p ' — —p 
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the expressions for Cx\i simplify to, 

Cii = 



(11) 



This remains true whether or not m — m' . If in addition it is also true that 
m — m' then (9) simplifies to 



u{p', A) ® u(p, - A) 



where 



-A-7o(i^ + m)(l - ATsf ) e^^"^ , (12) 



rj^ ^ rji - iXrj2. (13) 



These results also follow by using the relation [6] 

u(p, A) = -X-fou(-p, -A), 

in expression (4). 

Finally we remark that in the limit m — > equations (9) and (12) still apply 
with the proviso that 773 ^ p^ /m + 0{m/E). In the massless limit, where 
m\m —>■ 0, expression (9) simplifies to 
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u{p', X')u{p, A) = < 



1 



= Cva/^(1 + A75) for A' = -a, 

1' 

jC^i-xJ^l^f for A' = A. (14) 



[ 2^2^ 

It is easy to verify that (14) satisfies all of the identities valid for the inner 
product of two massless spinors. For example. 



and. 



m(p, -) m(p', +) = -(m(p, +) u{p\ -))*, 
\u{p, +)m(p', -)p = 2p-p'. 
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Equation (9). and its simplified forms (12) and (14), play a central role in our 
method. With them we can evaluate the trace in (3) and obtain an analytic 
expression for the amplitude in terms of dot products. 

Evaluating the trace, even at the amplitude level, may lead to very com- 
plicated expressions if the number of external gauge bosons connected to a 
single fermion line is large. Therefore, we develop an iterative method for 
obtaining this trace. The iterations can be easily performed either symboli- 
cally (using code hke REDUCE or FORM) or numerically. The basic idea is 
to use the fact that any product of an odd number of gamma matrices can 
be written in the form, 

W2---^2n+l=V'n + 754n- (15) 

Where the Vs and A's are obtained iteratively from, 

A': = 0. (16) 
This result follows after repeated use of the 7 identity^. 

For convenience we define the following expressions, 

Fx,yiBo,Vn,An) = uip,X)^^,(yn + lMuip',y). (17) 

The essence of the method is to write expressions of the form (1) in terms of 
F\,\'{yn, An) and Fx,a'(-Bo, Vn, An) defined above. The VnS and A^s are ob- 
tained by the simple iteration process (16). Finally, using (9) we express 
(17) in terms of four vector products. Since the general expressions for 

^In our conventions 75 = 17*^7^7^7^ and £0123 = 
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Fx^x'iVn, An) and Fx,a'(-^0) ^n) a-re common to every Feynman diagram 
we display their form below. 



where 



A",A 



Fx,x'{V,A) = 2AfJ2Cx',y'Mx",x 



A" 



Fx,x'{Bo,V,A) = 2AfJ2Cy,x"Nx",x 

A" 



m'p-V + mp'-V 

-Xlmm'-q^-A + rf-W{p,p', A, V)\ for A" = A 



mm'e^-A + e^-W{p,p\ A, V) 



(18) 



(19) 



for A" = -A, 



A",A 



and, 



' Bq ■ W{p, p', V, A) + mm'BQ ■ V 

+A Bq ■ W{rf, m'p + mp' , A, V) for A = A 

_ So • W{r]^, m'p + mp', A, V) for A" = -A, 



(20) 



W^'ix, y, z, w) = xf'y-z + zi^x-y - y'^x-z + ie'^^'^'x^ypWa- . (21) 



Of course these expressions take on much simpler forms when (12) and (14) 
are applicable. We will see an example of this simplification in the next 
section. We now have at hand all the tools necessary to express any object 
of the form (3) in terms of four vector products. In the next section we 
illustrate how these techniques are applied to compute the cross section for 
various processes. 



3 Applications 



As an illustration of our method we present the helicity amplitudes for two 
processes involving the top quark, namely, 
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Figure 1: Feynman diagrams for the process e^e — > tt 

and, 

77 — > tt. 

These processes wiU be important in the study of the phenomenology of the 
top quark at future electron-positron colliders. 

We can write the helicity amplitudes for both processes in terms of the spinor 
objects, 

A^(p,A;p',A') = «(p,A)7^t;(p',A') 

S^(p,A;p',A') = iZ(p,A)7^75^^(p',A') (22) 

We will work in the center of momentum system with the initial particles 
moving along the 2;-axis; we also take m = m'. In this frame we can then use 
(12) to obtain simple expressions for ^4^ and B^, 



-2mXg^iP' ifA' = A 
-2Er)^,ip) if A' = -A, 

(23) 

-2m if A' = A, 

-2Abl77^ ifA' = -A. 



Where rj^ip) are the polarization vectors defined in (13) corresponding the 
momentum p, mass m and helicity index A. 
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The e ^ tt Process: 



We now have the necessary ingredients to obtain expressions for the hehcity 
ampHtudcs. The Feynman diagrams for the first process are shown in figure 1. 
We take the 7^, Ztt, and Ze~e'^ couplings to be of the form'^, 



7ee 



67^ r^,,- = 7'^(a + 675). 
The amphtudes can then be written in terms of (23), 

T^(Ae, Ae-; A,, Af) = -^A^(Ae, Ae-)^ A^^iXt, Af) 



(24) 
(25) 



Tz{\e, Ag; At, At) — 



s — mt 



ac A^(Ae,Ae-)U'^(At,Aty 
+ ad [A^(Ae, Ag)^ S'^(At, Af)] (26) 
+ 6c [s^(Ae,Ag)tA''(At,At)" 
+ M [B^(Ae,Ag)tS''(At,Af; 



For simphcity we have suppressed the momentum indices, but the momentum 
dependence should be clear from the context. The electron and positron 
momenta are designated by Pe and pg, and the top, anti-top quark momenta 
are designated by p and p. In the CM frame the momenta take the form^. 



Ve = 


(E; 


0,0, E) 




Pe = 


{E; 


0,0, -E) 




P = 


(E; 


\p\ cos^,0, \p\ sinO) 






(E; 


-P), 


(27) 



*In the standard model the parameters qt, a, b, c, and d are given by 

4s2 - 1 



a = e ■ 



2e 



^ These are contravariant four vectors. 



9t = -y. 
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where E is the energy of the incoming electron in the center of mass frame. 
The electron mass is ignored in all of our expressions. Using these expres- 
sions for the momenta together with (23) and (26) it is simple to obtain an 
expression for the helicity amplitudes. We represent the total amplitude by, 
T — ^(Ae, Ae; At, Af), where the hehcity index for each particle is manifest. 
Then, 



r(Ae, -Ae; At, -At) = --^(AeAt + cos6') 

eqtPl + {a + bXe) {c + d Xt(3t) 

eqt(3l + {a + bXe)c 



T(Ae, — Ae; At, At) 



rrit Xt 
'E(3l 



(28) 



where (3z = yl — ml/AE"^, and /3t = Ip] /E. In obtaining the above expres- 
sion the following relations were useful. 



f]{-Xe,Pe) ■ ■n{Xt,p) = (cos 61 AeAt) 

p ■ fj{-Xe,Pe) = sine 

Pe-f]{-Xe,Pe) X ?7(At,p) = -i(At Ag COS 61) . 



(29) 



The 77 tt Process: 

As another example we consider the process 77 —>■ tt which will be present in 
the back-scattered laser mode of an e+e~ coUider. This experiment will have 
the possibility of polarizing both photons and of analyzing the polarization 
of the final state quarks. 

The relevant diagrams are those obtained via a t-quark exchanged in the 
i-channel (fig. 2). A straightforward application of (9) or (17) in the CM 
frame (27) leads to the following expression for the amplitude T, 



T{h,,h2,Xt,Xt)^^ 



T(]^){hi, /12, At, At) T^-!_-^{hi, /12, At, At) 



p-ki 



+ 



P ■ k2 



(30) 
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Figure 2: Feynman diagrams for the process 77 — > tt 



where A^, denote the top and anti-top hehcities and 'di^2 denote the photon 
polarization vectors, 



^1 = ^(1,^/^1,0); /2-;^(l,-i/i2,0), 
denoting by /ii,2 = ±1 the photon hehcities. We also defined 



Tg^{hi,h2,Xt,Xt)^2mt 
Tii^{h,,h2,Xt,-Xt)^2E 



Xf. 



\P\ 



(31) 



(32) 



In terms of the momenta in (27) we obtain, 
T{h,,h,Xt,Xt)^qt'^ ^ 



hi + Xt(5t (^AtAf 



T(/ii, -/ii. At, At) = gt ^2 2 

1 — Pt cos^ 



2 2 A A, sill i:^^''' 



if At = At, (33) 



'6 



{hi + Xtcos9) if At = —At- 



here f3t denotes the velocity of the top-quark, dt = \p\ /E, and 6 is the CM 
scattering angle. Note that as E 00, T(Xt, At, ±, =f) — as required by 
helicity conservation. Also in the non-relativistic limit the amplitude for two 
photons into two like-spin electrons vanishes, as it should. 



12 



4 Concluding Remarks 



We have presented a new method for computing hehcity amphtudes which is 
apphcable for processes involving both massless and massive fermions. The 
method allows one to obtain a covariant expression for the helicity ampli- 
tude without having to go through the process of squaring the amplitude. 
The method can be easily adopted to both numerical and symbolic compu- 
tations. We presented two examples which illustrate the simplicity of the 
method. Although for our two examples we worked in the center of mass 
frame the method is applicable in any frame and can be used in processes with 
many pairs of external massive fermions. In particular the methods presented 
here were used to compute the cross sections for gg —>■ ttH, gg ^ tt + bb, 
99 ~^ lltt, and many others (see for example references [14, 15]). For very 
complicated cases the combination of this method with a program for alge- 
braic manipulation (such as Form, Reduce, Macyma, Maple, or Mathemat- 
ica) should greatly simplify calculations. Even in those cases where one is 
only interested in the total cross section the use of the methods presented 
here may, in some cases, prove simpler to use than the traditional approach. 
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6 Appendix 

In this appendix we present a proof of the Michel-Bouchiat identity (4). We 
follow here the notation of reference [16]. 

The state of a particle can be characterized by giving the momentum p, 
mass m, total spin s, and its component of spin along a fixed direction. It is 
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common to take the 2;-axis or the direction of the momentum to specify this 
fixed direction. The latter case is refered to as the hehcity basis. In what 
follows we will use the helicity basis. The advantages of using the helicity 
basis are well known[l] and we will not dwell on them here. 

The hehcity operator is defined by, 

2W-h 

A = 

m 

where — ^e^^^^Pi^Jpa is the Pauli-Lubanski four vector; J^,y and are 
the angular momentum and energy-momentum operators respectively. The 
commutation relation for the W^s is, 

This relation follows from the commutation properties of Jij,i, and P^. The 
four vector n is any space-like normalized four vector with — —1. In the 
helicity basis n is in the direction of p. A conventional form for n is, 

\k\ h- 

m m\k\ 



Following Michel[17] we introduce an orthogonal set of tensors^, — {p^/m; 77^, 77^, 77^) 
such that, 

rf-p = 
rf-rf = -5'^ 

If we define, S°' = —W^U^/m, then it follows that, = S'-'q^^. For a particle 
of spin s, with rj & p • rj — 0, the operator, W • rj, has (2s + 1) eigenvalues: 
A = {— s, — s + 1, . . . , s — 1, s}. These eigenvalues are used to label states of 
given momenta and spin. 

From the commutation relation for the 1^'s and the properties of the 77's it 
is straightforward to show that the commutation relation for the 5"s obey an 
SU{2) algebra, 

[Si, Sj] — icijkSk, 
^The use of these tensors was first suggested to one of us by M.E. Peskin. 
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Furthermore, if we pick 77^ to have the form, 

^ \m m \k\) 

then up to a factor of |, the operator ^s, is just the hehcity operator. In 
the hehcity basis states are labeled by the eigenvalues of S^. The helicity 
projection operator is given by, 

= Hl±253), 

and as usual the lowering and raising operators of helicity are given by, 

S± — Si:k 182- 

We label the states for particles of spin s = i by uip, A) where, 

253m(p,A)= Am(p,A) 
^u{p,X)^ mu{p,X). (34) 

Then, 

u{p, ±) u(p, ±) = ii^^d^ + m), (35) 

where we are using a normalization such that, u{p, X)u{p, A) = 2m. 

In order to generalize the above expression to the case where the spinors have 
different helicity we simply use the raising and lowering operator for helicity, 

uiP,l')uip,±) = S^u{p,±)u{p,±) 
= St 2 ^ + 

where we have used, 5*:^ ^'3 = ±f 'S'q=. 

To obtain the Michel-Bouchiat form (4) we recall that. Si — —W • rji/m and 
use the definition of to obtain, 

2m ' 
15 



and, 

u{p,±)u{p,±) = l±^(^ + m), 

u{p,T)u{p,±) = ^,t!l^(j^ + m). 

Finally using the Pauli matrices, a*, these expressions can be written in the 
shorthand notation of (4). 
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